Docket Number STL11291 



EFFICIENT NOTCH COEFFICIENT COMPUTATION FOR A DISC 
DRIVE CONTROL SYSTEM USING FIXED POINT MATH 



Field Of The Invention 

The present invention relates generally to the use of notch filters in 
10 digital control systems of disc drives. In particular, the present invention 
provides a method and apparatus for efficiently calculating notch filter 
coefficients to compensate for the resonant frequencies of a particular disc 
drive. 

15 Background of the Invention 



Modern disc drive designs use digital feedback control systems to 
control the position of the head with respect to the tracks on the recording 
surface of the disc. Because disc drives and their electrical control systems, 

20 like all complex physical systems, are susceptible to resonance at particular 
frequencies, the feedback control system that is used to control the 
movement of the actuator arm generally employs one or more digital notch 
filters at some point in the feedback loop. These digital notch filters must 
be calibrated to attenuate the particular resonant frequencies of the 

25 particular disc drive in which the filters are used. In other words, a set of 
constant coefficients for each digital notch filter must be determined so that 
the filter attenuates an actual resonant frequency of the disc drive. 

Since the resonant frequencies of disc drive may vary from drive to 
drive, it is not sufficient to find a single set of coefficients for a given 

30 design. Rather, each disc drive must be individually calibrated. At present, 
this calibration process generally takes place at the time of media 
certification, which is when the manufacturer uses the assembled drive to 
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scan the recording media for defects and marks defective disc sectors as 
unusable. Manufacturers run a series of tests to determine the resonant 
frequencies of the drive. The disc drive, in turn, uses these resonant 
frequencies to look up filter coefficients in a table of coefficient values for 
various frequencies and programs its notch filters using the coefficients 
retrieved from the table. 

Table lookup, however, requires that all coefficient values be stored 
in memory. This makes it more complicated and costly to obtain precise 
filter coefficient values for a given design, since finer levels of precision 
require more firmware memory to store more filter coefficient values. 
What is needed, then, is a way to achieve precise filter calibration while 
minimizing memory usage. The present invention provides a solution to 
this and other problems, and offers other advantages over previous 
solutions. 
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Summary of the Invention 

The present invention provides a method and apparatus for 
efficiently calculating notch filter coefficients for a disc drive actuator arm 
control system. In a preferred embodiment, filter coefficients for a z- 
domain notch-filter transfer function are calculated in fixed-point 
arithmetic from polynomial interpolations of the non-linear functions that 
define the coefficients in terms of the notch frequency. These non-linear 
functions may be derived and interpolated a priori by applying the bilinear 
transform to an s-domain notch transfer function. Since, in a preferred 
embodiment, the z-domain transfer function can be expressed as a fraction, 
the numerator and denominator of the transfer function can be scaled so as 
to allow the coefficients to be expressed as integers, thus making it possible 
to calculate the filter coefficients from the aforementioned polynomial 
interpolations using fixed-point math. 
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Brief Description of the Drawings 

FIG. 1 is an exemplary perspective view of an exemplary disc drive; 

FIG. 2 is an exemplary top plan view of the printed circuit board of 
the exemplary disc drive of Figure 1; 

FIG. 3 is a block diagram of a control system for a disc drive in 
which the present invention may be utilized; 

FIG. 4 is a flowchart representation of a process of deriving 
interpolating polynomials for notch filter coefficients in accordance with a 
preferred embodiment of the present invention; and 

FIG. 5 is a flowchart representation of a process of deriving notch 
filter coefficients using fixed-point math in accordance with a preferred 
embodiment of the present invention. 
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Detailed Description 



Referring now to the drawings, and initially to FIG. 1, there is 
5 illustrated an example of a disc drive designated generally by reference 
numeral 20. Disc drive 20 includes a stack of storage discs 22a-d and a 
stack of read/ write heads 24a-h. Each of storage discs 22a-d is provided 
with a plurality of data tracks to store user data. As illustrated in FIG. 1, 
one head is provided for each surface of each of storage discs 22a-d such 

10 that data can be read from or written to the data tracks of all of the storage 
discs. The heads are coupled to pre-amplifier 31. It should be understood 
that disc drive 20 is merely representative of a disc drive system utilizing 
the present invention and that the present invention can be implemented in 
a disc drive system including more or less storage discs. 

15 Storage discs 22a-d are mounted for rotation by spindle motor 

arrangement 29, as is known in the art. Moreover, read/ write heads 24a-h 
are supported by respective actuator arms 28a-h for controlled positioning 
over preselected radii of storage discs 22a-d to enable the reading and 
writing of data from and to the data tracks. To that end, actuator arms 28a- 

20 h are rotatably mounted on pin 30 by voice coil motor 32 operable to 

controllably rotate actuator arms 28a-h radially across the disc surfaces. 

Each of read/ write heads 24a-h is mounted to a respective one of 
actuator arm 28a-h by a flexure element (not shown) and comprises a 
magnetic transducer 25 mounted to slider 26 having an air bearing surface 

25 (not shown), all in a known manner. As typically utilized in disc drive 

systems, sliders 26 cause magnetic transducers 25 of the read/ write heads 
24a-h to "fly" above the surfaces of the respective storage discs 22a-d for 
non-contact operation of the disc drive system, as discussed above. When 
not in use, voice coil motor 32 rotates actuator arms 28a-h during a contact 

30 stop operation, to position read/ write heads 24a-h over a respective one of 
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landing zones 58 or 60, where read/ write heads 24a-h come to rest on the 
storage disc surfaces. As should be understood, each of read/ write heads 
24a-h is at rest on a respective one of landing zones 58 or 60 at the 
commencement of a contact start operation. 

5 Printed circuit board (PCB) 34 is provided to mount control 

electronics for controlled operation of spindle motor 29 and voice coil 
motor 32. PCB 34 also includes read/ write channel circuitry coupled to 
read/ write heads 24a-h via pre-amplifier 31, to control the transfer of data 
to and from the data tracks of storage discs 22a-d. The manner for coupling 

10 PCB 34 to the various components of the disc drive is well known in the 

art, and includes connector 33 to couple the read/ write channel circuitry to 
pre-amplifier 31. 

Referring now to FIG. 2, there is illustrated in schematic form PCB 
34 and the electrical couplings between the control electronics on PCB 34 

15 and the components of the disc drive system described above. 

Microprocessor 35 is coupled to each of read/ write control 36, spindle 
motor control 38, actuator control 40, ROM 42 and RAM 43. In modern disc 
drive designs, the microprocessor can comprise a digital signal processor 
(DSP). Microprocessor 35 sends data to and receives data from storage 

20 discs 22a-d via read/write control 36 and read/ write heads 24a-h. 

Microprocessor 35 also operates according to instructions stored in 
ROM 42 to generate and transmit control signals to each of spindle motor 
control 38 and actuator control 40. 

Spindle motor control 38 is responsive to the control signals received 

25 from microprocessor 35 to generate and transmit a drive voltage to spindle 
motor 29 to cause storage discs 22a-d to rotate at an appropriate rotational 
velocity. 

Similarly, actuator control 40 is responsive to the control signals 
received from microprocessor 35 to generate and transmit a voltage to 
30 voice coil motor 32 to controllably rotate read/ write heads 24a-h, via 
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actuator arms 28a-h, to preselected radial positions over storage discs 22a- 
d. The magnitude and polarity of the voltage generated by actuator control 
40, as a function of the microprocessor control signals, determines the 
radial direction and radial speed of read/ write heads 24a-h. 

When data to be written or read from one of storage discs 22a-d are 
stored on a data track different from the current radial position of 
read/ write heads 24a-h, microprocessor 35 determines the current radial 
position of read/ write heads 24a-h and the radial position of the data track 
where read/ write heads 24a-h are to be relocated. Microprocessor 35 then 
implements a seek operation wherein the control signals generated by 
microprocessor 35 for actuator control 40 cause voice coil motor 32 to move 
read/ write heads 24a-h from the current data track to a destination data 
track at the desired radial position. 

When the actuator has moved read/ write heads 24a-h to the 
destination data track, a multiplexer (not shown) is used to couple 
read/ write heads 24a-h over the specific data track to be written or read, to 
read/ write control 36, as is generally known in the art. Read/ write control 
36 includes a read channel that, in accordance with modern disc drive 
design, comprises an electronic circuit that detects information represented 
by magnetic transitions recorded on the disc surface within the radial 
extent of the selected data track. As described above, each data track is 
divided into a number of data sectors. 

During a read operation, electrical signals transduced by the head 
from the magnetic transitions of the data sectors are input to the read 
channel of read/ write control 36 for processing via pre-amplifier 31. 
Random access memory (RAM) 43 can be used to buffer data read from or 
to be written to the data sectors of storage discs 22a-d via read/ write 
control 36. The buffered data can be transferred to or from a host computer 
utilizing the disc drive for data storage. 
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The present invention is directed toward a method and apparatus 
for generating notch filter coefficients for a disc drive actuator arm control 
system using fixed-point arithmetic. FIG. 3 is a diagram of a disc drive 
actuator arm control system with which a preferred embodiment of the 

5 invention may be employed. 

A digital command signal 300 for seeking or track-following is 
provided as input to summer 306, which combines command signal 300 
with feedback from plant 304, which in a preferred embodiment is voice 
coil motor 32 in FIG. 1. In a preferred embodiment, this feedback from 

10 plant 304 will be the position error signal (PES) of a particular disc drive 
head. A digital controller 301 generates a signal from the output of 
summer 306 that minimizes plant disturbance and maintains tracking 
requirements. Notch filter 302 is interposed between summer 306 and 
plant 304 and is used to filter out resonant frequencies from the control 

15 signal sent to plant 304, so as to prevent the control system in FIG. 3 from 
becoming unstable. In a preferred embodiment notch filter 302 may 
actually comprise a cascade of multiple filters configured to filter out a 
number of different resonant frequencies. 

A preferred embodiment of the present invention is implemented in 

20 the form of program code that is executed by a microprocessor within the 
disc drive. During the manufacturer's media certification process the 
resonant frequencies of the disc drive are ascertained by processes that are 
well-known to those skilled in the art. A preferred embodiment of the 
present invention uses these empirical resonant frequency measurements 

25 to derive coefficients for the actuator arm control system's notch filter. 

Although the calculations that are performed on the measured frequencies 
to obtain the coefficients are performed using fixed-point arithmetic, in 
order to best understand the principles of operation of a preferred 
embodiment of the present invention, it is helpful to delve into some 

30 mathematical analysis over the complex number system. 
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In analog signal processing, it is well-known and customary to 
express the transfer functions of system components as functions in the s- 
domain, i.e., as Laplace transforms. For example, the transfer function of a 
second-order notch filter is given as: 

5 

W + "8 (1) 

where coq is the notch frequency (expressed in radians per second) and Qi 
and Qi are "quality factors/' which are parameters that control the depth 
and width of the notch (stopband), respectively. The s-domain transfer 

10 function of a system is related to the frequency response of the system by 
the identity s=;co, where G) is a radian frequency and ;= ypi. 

Since the Laplace transform is a mathematical tool to performing 
analysis of continuous-time systems, the Laplace transform cannot be 
applied to digital (i.e., discrete-time) signal processing (at least not 

15 directly). Instead, the z-transform, which can be considered the discrete- 
time counterpart to the Laplace transform, is used. Since disc drives are 
generally digitally controlled, a preferred embodiment of the present 
invention uses fixed point arithmetic to calculate transfer function 
coefficients for a digital filter specified in terms of the discrete-time z- 

20 transform. 

In practice, it is customary in the design of digital filters to first start 
with an s-domain transfer function for the desired filtering function and 
transform the s-domain transfer function in to a z-domain transfer function. 
One commonly employed transformation technique for transforming an s- 
25 domain transfer function into a z-domain transfer function is known as the 
bilinear transformation (also referred to as Tustin's rule). The concepts of 
the s-domain, z-domain, and bilinear transformation are well-known in the 
digital signal processing art, so there is no need to provide an elaborate and 
detailed discussion of these topics. Nonetheless, a brief description of how 
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the bilinear transformation is performed, as well as some of the notation 
used in describing digital and analog signal processing concepts, may be 
helpful to an understanding of a preferred embodiment of the present 
invention. The reader seeking a more detailed treatment of the topic is 
directed to Ashok Ambardar, Analog and Digital Signal Processing, PWS 
Publishing, 1995, or any of the many introductory textbooks in the field of 
digital signal processing. 

A digital signal comprises a sequence of discrete values or samples, 
each representing a value of the signal at some point in time. A digital 
signal can be (and typically is) obtained from an analog signal by 
periodically sampling the magnitude of the analog signal over time to 
obtain the discrete digital signal values. The number of samples taken per 
second (i.e., the number of digital signal values representing a one-second 
interval of time) is referred to as the " sampling rate" of the digital signal, 

and is denoted as Sy. The sampling period t s =Sjr~^ measures the amount of 

time between successive samples. 

The concept of signal frequency is thus somewhat complicated for 
discrete-time signals, since there is both a sampling rate and whatever 
frequency components of the original sampled signal may exist. For 
example, a discrete-time representation of a sinusoidal function will have a 
signal frequency (i.e., the frequency of the sinusoid) as well as a sampling 
rate. The "digital frequency" of a signal is given by co / 1 $ , where 

co is the radian signal frequency. 

The bilinear transformation maps a function in the s-domain into the 
z-domain by performing the following substitution: 



In practice, however, an enhanced form of the bilinear 
transformation, the bilinear transformation with prewarping, is used to 
reduce frequency-related distortion. Prewarping is used to match the 
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response of the z-domain function at digital frequency £2rj> to the response 

of the s-domain function at radian signal frequency co. The bilinear 
transformation with prewarping utilizes the following substitution in place 
of that above: 

5 -> C . 

5 * + l (3) 

where C is defined as 

C = "A 

tan(|Q D ) ^ 

The bilinear transformation, therefore, allows a transfer function in 
the z-domain to be derived from an s-domain transfer function. It is well 

10 known in the art that transfer functions in the z-domain are readily 
implemented in digital circuitry. Indeed, z-domain transfer functions 
provide the basis for digital filters and digital signal processing in general. 

The present invention is concerned with the derivation of 
coefficients for z-domain transfer functions used for implementing digital 

15 notch filters in disc drives. Equation 1, above, provides the s-domain 

transfer function of a second-order notch filter. The bilinear transformation 
with prewarping can be used to transform this transfer function into a z- 
domain transfer function, as below: 

too - Ql — 2 * 

20 Multiplying the numerator and denominator by QiCMz+l) , 

multiplying out the parentheses, and collecting like terms yields 



a 2 z 2 j-_ai£+_ao 
[Z) ~ £> 2 z 2 +6^ + 60' (6) 



where 
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a 2 = Q X Q 2 C 2 + Q 2 wo + QiQ2^o, 
a x - 2Q 2 u; 0 + 2Q x Q 2 uSl - 2Q X Q 2 C 2 5 
a 0 = Q1Q2C 2 + QiQ 2 Wq + £2^0, 

fex = 2Q 1 ^ 0 + 2Q 1 Q 2 a; 0 2 -2g 1 Q 2 C 2 , 
60 = QiQ 2 C 2 + QiQ 2 ^ + Qi^ 0 - 

Thus, for a particular filter design, T(z) is the quotient of two polynomials 
in 2. The particular frequency response of the filter, then, is determined by 
the values of the coefficients of these polynomials. Thus, a primary task of 
a preferred embodiment of the present invention is to derive values for 
these coefficients using fixed-point math. 

There are two primary challenges that a preferred embodiment of 
the present invention overcomes in order to compute these coefficients 
using fixed-point math. First, since Q, C, and oq will generally have a 

significant fractional component, it is necessary to use scaling in order to 
allow the polynomial coefficients to be expressed as integers without losing 
a great deal of accuracy. Since T(z) is a fraction, Equation 6 can be scaled 
by simply multiplying the numerator and denominator by a constant K, 
where the choice of K will depend upon the general frequency range of 
interest (i.e., it is a design choice that can be made in advance). 

The second challenge that a preferred embodiment of the present 
invention overcomes relates to the fact that in Equations 4 and 6, C is a 
nonlinear function of Qp, the digital frequency, since C contains the 

trigonometric tangent function. Computing a tangent function requires an 
ability to compute fractions or decimals. Thus, it is not possible to apply 
Equation 6 directly using fixed-point math. 

A preferred embodiment of the present invention, instead, uses 
polynomial interpolation to derive a value for each coefficient in Equation 
6. There are a number of different techniques known in the art for deriving 
a polynomial interpolation of a function. One particular technique that is 
well-known in the art is known as Lagrange polynomial interpolation. In 
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Lagrange polynomial interpolation, to approximate a function using an 
nth-degree polynomial, one first evaluates the original function, say f(x), at 
n+1 points (xQ,xy.,.,x n ) in order to obtain a series of data points, 

f( x 0)'f( x l)'-- 'f( x n)- The "til-degree polynomial interpolation is given by 

n 

Pn(x) = J2 l i( X )f( X i)' 

*=0 (7) 

where the l(s are given by 

u(x)= n (^^). 

,4j=o V^i- xjJ (8) 
One of ordinary skill in the art will recognize that the expressions 
given above for Lagrange polynomial interpolation can be expanded 
symbolically so that the expanded form of p n can be determined by simply 

arithmetic evaluation in a computer. This, and other methods of 
polynomial interpolation, are well-known in the art as being readily 
performed by a computer. The interested reader may consult any suitable 
undergraduate textbook in numerical analysis for further details regarding 
polynomial interpolation and other function-approximation techniques. 

In the context of a preferred embodiment of the present invention, a 
polynomial interpolation of each of each of the coefficient expressions in 
the filter's transfer function is made. Specifically, each of the coefficient 
expressions in the transfer function is expressed as a function of the notch 
frequency g)q, and an interpolating polynomial in odq is obtained. 

For example, consider the expression given in Equation 6 for the 

coefficient of in the numerator of the transfer function. We will then 
treat this expression as a function /of the desired radian notch frequency go, 
as below: 



14 



Docket Number STL11291 



f(u>) - Q1Q2C 2 + Q 2 UJ + (9) 
If we make the a priori design decisions to match the digital and 
analog transfer functions at the same radian notch frequency, use a 
sampling frequency of 13.32 kHz, and to use quality factors of Qi= Q2=l, 
we can substitute Equation 4 into Equation 9 to obtain 

/H = ( tan(6660a,) ) + " + ' (10 ) 
Using floating-point math, we can evaluate /(co) at n+1 points to 
obtain /(co Q),f((o \)>-A' s> «)• In a preferred embodiment, n=2. To continue 
with the example, let us suppose that we evaluate Equation 10 at co 
0 =27r-440, co 2=2^-880, and co 2 =2^-1760. We then obtain /(co 0 )«7514402.33, 
/(co 2) «30446290.00, and/(co 2 ) «1221 68311. 57. Polynomial interpolation can 

then be performed on these values using either floating point math or by 
using fixed point math on rounded or truncated versions of the values. For 
example, Lagrange polynomial interpolation would proceed as below: 

/u-cj 1 \/u}-oj 2 \_ f u) - 27T • 880 \ / u> - 2it ■ 1760 N 
l °( u ) ~ [coq-ujJ \u, 0 -u 2 ) ~ V2tt • 440 - 2tt • 880 J \2ir ■ 440 - 2tt - 17607 

(11) 



uj - 2?r • 440 \ / oj - 2tt • 1760 



/ l \ - ( UJ - UJ o \ ( u-u) 2 \ _ / u - 2tt • 440 \ / 
lta,j ~ Ui - ujJ Ui - wj V27T • 880 - 2tt • UOJ V 



2tt • 880 - 2?r • 1760 
(12) 



u - 27r • 440 \ / uj - 2tt • 880 



, / x _ / u; - u; 0 \ / cj - cji \ = / u)-2tt- 440 \ / 

1 2 (oj) - ^ _ ^ _ w j ^ 176Q _ 27r . 440 A2tt • 



1760 - 2?r • 880 
(13) 



P(w) = 13 k(w)f{wi) (14) 

i=0 

= Zo(w)/(wb) + ii(w)/(wi) + fa(w)/W (15) 

ss cj 2 + 30437996a; - 84149059540. ( 16 ) 
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Equation 16, above, shows that the interpolating polynomial p(co) in 
our example can be evaluated using fixed-point math, since the only 
operations that are needed in order to evaluate p(co) are addition, 
5 subtraction and multiplication over the integers. Moreover, even if it is 
necessary to scale the z-domain transfer function from which p(co) is 
derived (i.e., by multiplying the numerator and denominator of the transfer 
function by a constant k), p(co) can be easily scaled as well, simply by 
multiplying each term of /?(co) by k. Moreover, one of ordinary skill in the 

10 art will recognize that a scaled form of co, for example co '=fc co (k being the 

scale factor), may also be used by simply substituting co'//c for co in p(co). 

Thus, using the above-described techniques, each coefficient of a z- 
domain notch filter can be expressed as a polynomial over the integers that 
is a function of the desired radian notch frequency co. Hence, appropriate 

15 coefficients for a notch filter having a particular notch frequency may be 
calculated by a computer by supplying the notch frequency as input and 
evaluating the appropriate polynomial approximation for each coefficient 
in the z-domain transfer function of the filter. In a preferred embodiment, 
the interpolating polynomials for each coefficient are derived a priori using 

20 floating-point math, and the disc drive is programmed to evaluate these 
polynomials using fixed-point math. 

FIG. 4 is a flowchart representation of a process of deriving 
interpolating polynomials for use in calculating notch filter coefficients in 
accordance with a preferred embodiment of the present invention. The 

25 bilinear transformation is used to take an s-domain transfer function and 
convert it into an equivalent z-domain transfer function (block 400). The z- 
domain transfer function so obtained is scaled so as to make the transfer 
function coefficients expressable as integers without significant loss of 
precision (block 402). Finally, polynomial interpolation is performed on 

30 each of the expressions defining the coefficients in the z-domain transfer 
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function, so as to obtain interpolating polynomials that approximate the 
desired transfer function coefficients as functions of the filter's notch 
frequency (block 404). 

FIG. 5 is a flowchart representation of a process of assigning 

5 appropriate notch filter coefficients to a disc drive control system notch 
filter in accordance with a preferred embodiment of the present invention. 
A resonant frequency of the disc drive control system is first determined 
(block 500); in the case of several resonant frequencies, each is determined 
and the subsequent steps are performed with respect to each frequency, so 

10 that a separate notch filter is derived for each notch frequency. Next, the 
approximating polynomials for each of the transfer function coefficients are 
evaluated at the resonant frequency (i.e., the desired notch frequency for 
the filter) (block 502). Finally, the results of evaluating the approximating 
polynomials are used as the filter coefficients to program the digital notch 

15 filter (block 504). 

It is important to note that while the present invention has been 
described in the context of a fully functioning data processing system, 
those of ordinary skill in the art will appreciate that the processes of the 
present invention are capable of being distributed in the form of a 

20 computer readable medium of instructions or other functional descriptive 
material and in a variety of other forms and that the present invention is 
equally applicable regardless of the particular type of signal bearing media 
actually used to carry out the distribution. Examples of computer readable 
media include recordable-type media, such as a floppy disk, a hard disk 

25 drive, a RAM, CD-ROMs, DVD-ROMs, and transmission-type media, such 
as digital and analog communications links, wired or wireless 
communications links using transmission forms, such as, for example, 
radio frequency and light wave transmissions. The computer readable 
media may take the form of coded formats that are decoded for actual use 

30 in a particular data processing system. Functional descriptive material is 
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information that imparts functionality to a machine. Functional descriptive 
material includes, but is not limited to, computer programs, instructions, 
rules, facts, definitions of computable functions, objects, and data 
structures. 

5 The description of the present invention has been presented for 

purposes of illustration and description, and is not intended to be 
exhaustive or limited to the invention in the form disclosed. Many 
modifications and variations will be apparent to those of ordinary skill in 
the art. The embodiment was chosen and described in order to best 

10 explain the principles of the invention, the practical application, and to 
enable others of ordinary skill in the art to understand the invention for 
various embodiments with various modifications as are suited to the 
particular use contemplated. 



